Morphology of the nonspherically decaying radiation generated 
by a rotating superluminal source: reply to comment 
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Abstract 

The fact that the formula used by Hannay in the preceding Comment is "from a standard text 
on electrodynamics" neither warrants that it is universally applicable, nor that it is unequivocally 
correct. We have explicitly shown [J. Opt. Soc. Am. A 25, 543 (2008)] that, since it does not 
include the boundary contribution toward the value of the field, the formula in question is not 
applicable when the source is extended and has a distribution pattern that rotates faster than light 
in vacuo. The neglected boundary term in the retarded solution to the wave equation governing 
the electromagnetic field forms the basis of diffraction theory. If this term were identically zero, 
for the reasons given by Hannay, the diffraction of electromagnetic waves through apertures on a 
surface enclosing a source would have been impossible. 



I. INTRODUCTION 

The argument presented by Hannay in Ref. Q is based on an incorrect solution of Maxwell's 
equations. We have explicitly shown that the retarded solution to the wave equation 

V 2 B-^ = - T Vxj. (1) 

governing the magnetic field B is not always given by 

B(x P) f P ) = ± / d 3 J^L (2) 
c J |x P - x| 

as assumed by Hannay [l!]; an exception is the solution of Maxwell's equations describing 
the emission from a polarization current density j whose distribution pattern rotates su- 
perluminally (i.e., faster than light in vacuo). [Here, (xp,£p) and (x, t) are the space-time 
coordinates of the observation point and the source points, respectively, c is the speed of 
light in vacuo, and the square brackets denote the retarded value of V x j.] The emission 
from such a source consists of a collection of narrowing subbeams for which the absolute 
value of the gradient of the radiation field B increases (as Rp 7 ^ 2 ) with the distance Rp from 
its source, rather than decreasing as predicted by Eq. (jSJ) [2j. The inadequacy of Eq. (j2J), in 
this case, lies in the neglect of the boundary contribution toward the value of B. 

We first describe how the boundary contribution to the retarded solution for the potential 
can always be made equal to zero, irrespective of the source motion. In the Lorenz gauge, 
the electromagnetic fields 

E = -VpA° -dA/d(ctp), B = V P xA, (3) 

are given by a four-potential N 1 that satisfies the wave equation 

1 8 2 A^ Air 

^ 2 ^-- 2 ^ = --f, A* = 0,-,3, (4) 
cr ot c 

where A / c and j / c are the electric potential and charge density and A M and j M for \i — 1,2,3 
are the Cartesian components of the magnetic potential A and the current density j jsj. The 
solution to the initial-boundary value problem for Eq. is given by 

A"(xp, tp) =- [ P dt [ d 3 x fG + — f V dt [ dS • (GVA M - A^VG) 

c Jo JV 47r Jo Jy, /c\ 

1 d 3 x[A^— G » 



in which G is the Green's function and £ is the surface enclosing the volume V (see, e.g., 
page 893 of Ref.|j). 

The potential arising from a general time- dependent localized source in unbounded space 
decays as -Rp" 1 when Rp = |xp| — > oo, so that for an arbitrary free-space potential the 
second term in Eq. §5§ would be of the same order of magnitude (~ -Rp -1 ) as the first 
term, in the limit that the boundary £ tends to infinity. However, even potentials that 
satisfy the Lorenz condition V ■ A + c~ 2 dA° jdt = are arbitrary to within a solution of the 
homogeneous wave equation: the gauge transformation 

A^A + VA, A -> A -dA/dt (6) 

preserves the Lorenz condition if V 2 A — c~ 2 d 2 A/dt 2 = (see Ref. 0). One can always use 
this gauge freedom in the choice of the potential to render the boundary contribution (the 
second term) in Eq. (j^J) equal to zero, since this term, too, satisfies the homogeneous wave 
equation. Under the null initial conditions A^\ t= Q = (cL4 M /dt) t =o = assumed in this note, 
the contribution from the third term in Eq. (jSJ) is identically zero. 

In the absence of boundaries, the retarded Green's function has the form 

x 5(tp - t - R/c) 
G(x, t; x P , tp) = lF (7) 

where 5 is the Dirac delta function and R is the magnitude of the separation R = xp — x 
between the observation point xp and the source point x. Irrespective of whether the 
radiation decays spherically (as in the case of a conventional source) or nonspherically (as 
would apply for a rotating superluminal source) [2|, therefore, the potential A^ due to a 
localized source distribution that is switched on at t = in an unbounded space, can be 
calculated from the first term in Eq. (JHJ): 

A M (xp, t P ) = c- 1 J d 3 xdt j"(x, t)S(t P - t - R/c)/R, (8) 

i.e., from the classical expression for the retarded potential. Whatever the Green's function 
for the problem may be in the presence of boundaries, it would approach that in Eq. (j7j) in 
the limit where the boundaries tend to infinity, so that one can also use this potential to 
calculate the field on a boundary that lies at large distances from the source. 

We now turn to the case of the field and show that an analogous assumption about 
the boundary contribution may not be made. Consider the wave equation ([!]) governing the 



magnetic field; Equation ([I]) may be obtained by simply taking the curl of the wave equation 
for the vector potential [Eq. PJ for [A — 1, 2, 3]. We write the solution to the initial-boundary 
value problem for Eq. ([T]), in analogy with Eq. §5§, as 

£? fc (x P , t P ) =- [ P dt [ d 3 x (V x j) k G + f P dt f dS • (GVB k - B k VG) 

c Jo Jv 47r Jo </£ (q\ 



47T C 2 Jy V dt dt ) t=0 

where k = 1, 2, 3 designate the Cartesian components of B and V x j. 

Here, we no longer have the freedom, offered by a gauge transformation in the case of Eq. 
flSJ), to make the boundary term zero. Nor does this term always decay faster than the source 
term, so that it could be neglected for a boundary that tends to infinity, as is commonly 
assumed in textbooks (e.g., page 246 of Ref. y|) and the published literature 
The boundary contribution to the retarded solution of the wave equation governing the 
field [the second term on the right-hand side of Eq. (j2j)] entails a surface integral over the 
boundary values of both the field and its gradient. For a rotating superluminal source, 
where the gradient of the field increases as Rp 7 ^ 2 over a solid angle that decreases as -Rp -4 , 
this boundary contribution turns out to be proportional to Rp~ 1 ^ 2 (see Ref. |2|). Not only is 
this not negligible relative to the contribution from the source term [the first term on the 
right-hand side of Eq. (jUJ)], but the boundary term constitutes the dominant contribution 
toward the value of the radiation field in this case. 

If one ignores the boundary term in the retarded solution of the wave equation governing 
the field (as Hannay does , one obtains a different result, in the superluminal 



ne sup 

0,Q, 



regime, from that obtained by calculating the field via the retarded potential [9|,ll0|,lll|. This 
apparent contradiction stems solely from having ignored a term in the solution to the wave 
equation that is by a factor of the order of Rp 1 ^ 2 greater than the term that is normally kept 
in this solution. The contradiction disappears once the neglected term is taken into account: 
the solutions to the wave equations governing both the potential and the field predict that 
the field of a rotating superluminal source decays as Rp^ 1 ^ 2 as Rp tends to infinity, a result 
that has now been demonstrated experimentally jl2j |. 

For a volume V in which there are no sources, Eq. (Q reduces to 

fl fc (x P ,tp) = ^ / P dt / dS • (GVB k - B k VG), (10) 

J J s 
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under the null initial conditions B k \ t= o = (dB k /dt) t =o = 0. As in the customary geometry 
for diffraction, the closed surface S can consist of two disjoint closed surfaces, £ inner and 
Souter (e.g., two concentric spheres), both of which enclose the source (see Fig. 10.7 of Ref. []). 
If the observation point does not lie in the region between S inner and S outer , i.e., lies outside 
the closed surface S, then the composite surface integral in Eq. (TTU|) vanishes: 



f F dt j dS • (GVB k - B k VG) = f P dt( f + I \ dS • (GVB k - B k VG) 

JO JT. Jo VE 

inner ^outcr / 
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fsee Ref. 



4|). But under no circumstances would the integrals over Sinner or S outC r vanish 
individually if these surfaces enclose a source, i.e., if there is a nonzero field inside Shiner- Nor 
does the fact that the values of these integrals are unchanged by deformations of S inn er and 
S ou ter have any bearing on whether they are nonvanishing or not. Equation ffTOl) forms the 
basis of diffraction theory [3]. If the surface integrals over S inn er and S outor were to vanish 
individually, as claimed by Hannay 1|, the diffraction of electromagnetic waves through 
apertures on a surface enclosing a source would have been impossible (see Sec. 10.5 of 
Ref. J. 

Finally, we must stress that what one obtains by including the boundary term in the 
retarded solution to the wave equation governing the field is merely a mathematical identity; 
it is not a solution that could be used to calculate the field arising from a given source 
distribution in free space. Unless its boundary term happens to be small enough relative to 
its source term to be neglected, a condition that cannot be known a priori, the solution in 
question would require that one prescribe the field in the radiation zone {i.e., what one is 
seeking) as a boundary condition. Thus, the role played by the classical expression for the 
retarded potential in radiation theory is much more fundamental than that played by the 
corresponding retarded solution of the wave equation governing the field. The only way to 
calculate the free-space radiation field of an accelerated superluminal source is to calculate 
the retarded potential and differentiate the resulting expression to find the field (see also 



Ref. 
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